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Abstract: The renormalization group (RG) approach to cosmology is an efficient method 
to study the possible evolution of the cosmological parameters from the point of view of 
quantum field theory (QFT) in curved space-time. In this work we continue our previous 
investigations of the RG method based on potential low-energy effects induced from physics 



CO 



at very high energy scales Mx ^ Mp. In the present instance we assume that both the 



Newton constant, G, and the cosmological term, A, can be functions of a scale parameter \x. 
It turns out that G(fi) evolves according to a logarithmic law which may lead to asymptotic 
freedom of gravity, similar to the gauge coupling in QCD. At the same time A(/x) evolves 
quadratically with /u. We study the consistency and cosmological consequences of these 
laws when [i ~ H. Furthermore, we propose to extend this method to the astrophysical 
domain after identifying the local RG scale at the galactic level. It turns out that Kepler's 
third law of celestial mechanics receives quantum corrections that may help to explain the 
flat rotation curves of the galaxies without introducing the dark matter hypothesis. The 
origin of these effects (cosmological and astrophysical) could be linked, in our framework, 
to physics at M x ~ 10 16 ~ 17 GeV. 
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1. Introduction 



The recent developments in physical cosmology have provided an astonishingly accurate 
picture of our universe, within the canons of the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) paradigm From these achievements emerged what has been called the cos- 
mological concordance model, characterized by an essentially zero value of the spatial cur- 
vature parameter and a non-vanishing positive value of the cosmological constant (CC) 
term, A, in Einstein's equations - or in general of some form of dark energy (DE) that 
cannot be attributed to any form of matter and radiation. This new standard cosmological 
model supersedes and refutes (by many standard deviations) the old Einstein-de Sitter (or 
critical-density) cosmological model, which has been the standard cosmological model until 
recently. The latter is also a spatially flat model, but it is characterized by a zero value of 
the cosmological constant. Evidence for the values of the cosmological parameters within 
the new concordance model comes both from tracing the rate of expansion of the universe 
with high-z Type la supernovae (SNe la) experiments ||, ||| and from the precise mea- 
surement of the anisotropies in the cosmic microwave background (CMB) radiation ||, |||] . 
These experiments indicate that the value of the cosmological constant is around 1 

A = n° A p° c ~ Qhl x W~ 47 GeV 4 , (1.1) 

where ~ 70%. Here p® ~ (3.0y/ho x 10~ 12 GeV) 4 is the present value of the critical 
density and Hq = 100 /io km/s/Mpc, with ho = 0.72 db 0.08, is the most precisely known 
value of Hubble's constant ||. In the context of the Standard Model (SM) of electroweak 
interactions, the physical (measured) value ( |1.1| ) should be the sum of the vacuum CC term 
in Einstein's equations, A mc , and the induced contribution from the vacuum expectation 
value of the Higgs effective potential, A in d = (V c g), namely A = A vac + A inc [ . Other 
contributions should be added within the SM, like those from QCD, but they are much 
smaller. It is only the combined parameter A that makes physical sense, whereas both A vac 
and Ai n d remain individually unobservable- see e.g. m for a more detailed explanation. 
From the current LEP numerical bound on the Higgs boson mass, Mh > 114.4 GeV (at the 
95% C.L.) f], one finds \A ind \ = (l/8)M|i? 2 > 1.0 x 10 8 GeV 4 , where v ~ 246 GeV is the 
vacuum expectation value of the Higgs field. Clearly, | Ai n< i\ is 55 orders of magnitude larger 
than the observed CC value ([O]). Such discrepancy, the so-called "old" CC problem || [l(|, 
manifests itself in the necessity of enforcing an unnaturally exact fine tuning of the original 
cosmological term A vac in the vacuum action that has to cancel the induced counterpart 
Ai n d within a precision (in the SM) of one part in 10 55 . But the CC problem seems to be 
quite more general and it seems to plague all physical theories (QFT and string theory) 
making use of the concept of vacuum energy |], [l0|, p]] . Moreover, all attempts to deduce 
the small value of the cosmological constant from a sound theoretical idea involve some 
explicit or implicit form of severe fine-tuning among the parameters of the model. This 
became clear already from the first attempts to treat the dark energy component as a 

1 The original CC term, A, in Einstein's equations is related to our A by A = 8ir G A, where G is Newton's 
constant. Notice that A, in our notation, has dimensions of energy density, and to emphasize this fact it 
will sometimes be denoted K = p\. 
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dynamical scalar field 12, 13, 14]. More recently this approach took the popular form of 



a "quintessence" field slow-rolling down its potential [15| and has adopted many different 



forms [16, 17]. The main motivation for the quintessence models is that they could provide 
a variable dark energy. Although there is no compelling reason for a time evolving DE at 
the moment - and it is difficult to get an experimental handle on it [^]- it is even more 
difficult to accept (at least from the theoretical point of view) that the DE can be described 
by a static CC term in Einstein's equations with the small value fll.ip throughout most of 
the history of the universe, namely after the electroweak phase transition took place. 

An alternative possibility is to think of the DE as a time-evolving CC. This has been put 



forward from a purely phenomenological point of view in many places of the literature [19 



2C|. However, a treatment of an evolving CC term from the more fundamental point of 
view of the renormalization group in QFT, with a direct link to Particle Physics, has been 
suggested only recently |2l], 0, |22|, The possible experimental consequences of this idea 
in the light of high-z SNe la data have also been tested in great detail in |24| (see also the 
framework and analysis of |25| , |26||). In the renormalization group method one takes the 
point of view that A should not be constant, because the quantum effects may shift away 
the prescribed value, even if the latter is assumed to be zero|2l|. In this approach the 
equation of state for the DE is the one for the "true" cosmological term p\ = —p\. Hence 
it is fundamentally different from all kinds of quintessence-like proposals, where the CC is 
assumed to vanish precisely and the role of the DE is played by some new entity with a 
nontrivial equation of state. In the RG approach the CC becomes time-evolving, or redshift 
dependent: A = A(z). Although we do not have a quantum theory of gravity where the 
running of the gravitational and cosmological constants could ultimately be substantiated, 
a semiclassical description within the well established formalism of QFT in curved space- 



time (see e.g. [27, [28|]) should be a good starting point. From the RG point of view, A 
becomes a scaling parameter A(/i) whose value should be sensitive to the entire energy 
history of the universe - in a manner not essentially different to, say, the electromagnetic 
coupling constant in QED or the gauge coupling constant in QCD. 



While the RG method to tackle the CC problem has been used earlier [29, [}(], gjj , 
the decoupling effects of the massive particles at low energies may change significantly the 
structure of the renormalization group equations (RGE) with important phenomenological 
consequences. This feature has been pointed out recently by several authors from various 
interesting points of view p2|, [jj. However, it is not easy to achieve a RG cosmological 



model where A runs smoothly without fine tuning at the present epoch. In Ref. [23, 32] 



a successful attempt in this direction has been made, which is based on possible quantum 
effects near a high mass scale Mx, not far away from the Planck scale Mp, that are trans- 
ferred (via soft decoupling) to the low energy physics. At the same time, the approximate 
coincidence of the observed CC and the matter density, ~ 09, (i.e. the "new" CC 
problem, or "time coincidence problem" j|, [R]]) can be alleviated in this framework p2|| . 
Thus the RG cosmologies constitute a healthy alternative to the ad hoc quintessence models 
and are worth exploring 2 . In the present paper we generalize the RG approach developed 



2 See |33| for a summarized presentation of this RG framework and an exposition of the basic results. 
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in jj], 22, 32 1 by including the possibility to have a running Newton's constant as well 3 . 
A fundamental question that we address here is the possibility to maintain the phenomeno- 
logically successful quadratic evolution law 5 A ~ H 2 (emerging from an RGE of the form 
dA/dlnH oc H 2 ) as proposed in |7], 23, 32], also in a G-running framework. This law is 
based on the assumption that the physical RG scale is defined by fi = H, the Hubble 
parameter. We will see that this Ansatz is consistent with the systematic RG scale-setting 
procedure proposed recently in In addition, the 5 A ~ H 2 law is the only one among 
those that are consistent with this procedure which could have measurable phenomenolog- 
ical consequences. Therefore, it is of foremost importance to explore its phenomenological 
viability, in spite of the insurmountable technical difficulties encountered at present for an 
explicit calculation in QFT on a curved background |3?]] . In this sense it is remarkable that 
the quadratic law 5 A ~ H 2 can also be motivated from very general arguments in QFT [38] 
which may go beyond the RG. No less appealing is that this quadratic law can further be 
claimed, at least heuristically, on the grounds of the so-called Holographic Principle (HP) 
applied to effective field theories In fact, by linking IR and UV scales through entropy 
bounds, the HP may provide essential ingredients to relate the microscopic world of Par- 
ticle Physics and the large scale structure of the Universe 4 . A main result of this paper 
is that the adoption of that quadratic law will lead us to an alternative cosmological RG 
framework to that presented in Ref. (2^, [32|, keeping though all the nice phenomenological 
features of the latter, however with the distinctive property of having a (logarithmically) 
running gravitational constant. This last property endows the new cosmological model 
with potentially far-reaching phenomenological implications both at the cosmological and 
astrophysical level which will be explored in some detail in this paper. 

The structure of the paper is as follows. In the next section we introduce a few general 
features of the cosmological models based on variable G and A. In Section [|| we provide 
some general discussion of the RG cosmological models. The consistency of the RG scale 
setting in cosmology is discussed in Section |4|. In Section [B| we focus on the cosmological 
model presented here, characterized by a logarithmic running of G and a quadratic running 
of A. The numerical analysis of this model is presented in Section || where we also show the 
predicted deviations from the standard FLRW expectations. In Section [?] we narrow down 
the applicability of this RG cosmological model to the astrophysical domain by considering 
the possibility to attribute the origin of the flat rotation curves of the galaxies to a RG 
quantum correction to Kepler's third law of celestial mechanics. We close with some final 
discussion in the last section, where we draw our conclusions. 



2. Variable G and A cosmologies 

The cosmological constant contribution to the curvature of space-time is represented by 
the Xg^v term on the l.h.s. of Einstein's equations. The latter can be absorbed on the 

Recently some applications of the RG method to inhomogeneous cosmologies have also been proposed [0 . 

3 Phenomenological models with variable G and A, without reference to the RG, have been studied in 
the literature since long ago ^(J . See also the more recent papers J3f| . 

4 See e.g. Mf for a general review of the HP and its applications. 
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r.h.s. of these equations 

R„u-^g f ,uR = 87TGf llu , (2.1) 

where the modified is given by T^ v = T^ u + g^ v A. Here T^ u is the ordinary energy- 
momentum tensor associated to isotropic matter and radiation, and the new CC term 
A = A/87T G = pa represents the vacuum energy density. Modeling the expanding universe 
as a perfect fluid with velocity 4- vector field U^, we have 

V = -pg^ + {P + P)U^U U , (2.2) 

where p is the proper isotropic pressure and p is the proper energy density of matter- 



radiation. Clearly the modified T^ v defined above takes the same form as (2^2) with 
p — > p = p + A and p — > p = p — A. With this generalized energy-momentum tensor, and 
in the FLRW metric (k = for flat, k = ±1 for spatially curved, universes) 

ds 2 = dt 2 _ a 2f t) ( + r 2 dQ 2 + r 2 gin 2 Q ^ fr y 

\ 1 — k r z ) 
the gravitational field equations boil down to Friedmann's equation 
Tr2 fa\ 2 8irG k 8vrG, k 



and the dynamical field equation for the scale factor: 



-— G(p + 3p)a = -— G(p + 3p- 2A)a. (2.5) 
3 3 



The measured CC is the parameter A entering equations (2.4) and (|2.5|), which are used to 



fit the cosmological data ||, [J. Let us next contemplate the possibility that G = G(t) and 
A = A(t) can be both functions of the cosmic time. This is allowed by the Cosmological 
Principle embodied in the FLRW metric ([2.3|). Then the Bianchi identities imply that 
y ^ I G T^u I = 0, and with the help of the FLRW metric a straightforward calculation 



yields 

j [G{A + p)\ +3GH (p + p) =0. (2.6) 

This equation just reflects the purely geometric Bianchi identity satisfied by the tensor on 
the l.h.s. of Einstein's equations ( [2.1| ) in terms of the physical quantities (energy densities, 
pressure and gravitational coupling) on its r.h.s. It is easy to check that ( |2.6| ) constitutes 
a first integral of the dynamical system (^4|) and (^1]). However, this general constraint 
on the physical quantities is too complicated in practice, and it is difficult to understand 
it on physical grounds because it does not generally lead to a local energy-momentum 
conservation law. Let us consider some simpler cases. The simplest possible case is when 
both G and A are constant. When G is constant, the identity above implies that A is also 
a constant, if and only if the ordinary energy-momentum tensor is individually conserved 
(V M ^ = 0), i.e. 

p + 3H (p + p) =0. (2.7) 
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However, a first non-trivial situation appears when A = A(i) but G is still constant. Then 
a first integral of the system ( |2.4| ) and ( |2.5| ) is given by 



A + p + 3#(p + p) = 0. (2.8) 

This result is simpler than (^1]), and has the advantage that it can still be interpreted 
as a physical conservation law. The conserved (total) energy density here is the sum of 
the matter-radiation energy density, p = pu + PR, plus the dark energy density embodied 
in pa = A. Of course pm includes the contribution from the purported cold dark matter 
(CDM), which should be dominant. Alternatively, using the metric ( |2.3| ) one may derive 
(|2.Sj ) directly from y M = 0, as this is obviously a particular case of (|2.6|) . This possibility 
shows that a time- variable CC cosmology may exist such that transfer of energy occurs from 
matter-radiation into vacuum energy, and vice versa. From the point of view of general 
covariance there is no problem because the latter is tied to the conservation of the full 
energy-momentum tensor, which in this case is T^ u rather than T^ u . We remark that while 
the decay of the vacuum energy into ordinary matter and radiation could be problematic 



for the observed CMB, the decay into CDM could be fine [16, 41]. The detailed analysis 
of a time- variable CC cosmology in the light of high-z SNe la data carried out in Ref. [24] 
was precisely based on this framework. 

In this paper we make a step forward by entertaining the possibility that G can also 
be a function of the cosmic time, G = G(t), together with A = A(i). Although we have 
seen that the most general case should fulfill the general constraint fl2,6| ), we still want to 
keep the usual matter-radiation conservation law fl2.7| ) as this represents the canonical (i.e. 
the simplest) cosmological model with variable G and A possessing a physical conservation 
law. In this model there is no transfer of energy between matter-radiation and dark en- 
ergy. However, the time evolution of A is still possible at the expense of a time varying 
gravitational constant. Indeed, it is easy to see that in this case the following (hereafter 
called canonical) differential constraint must be satisfied: 

(p + A)G + GA = 0. (2.9) 



Again this is a particular case of fl2.6| ). The solution of a generic cosmological model of 
this kind with variable G and A is contained in part in the coupled system of differential 
equations fl2.4|) , ( |2.7D and ( |2.9[ ), together with the equation of state p = p{p) for matter 
and radiation. However, still another equation is needed to completely solve this cosmolog- 
ical model in terms of the fundamental set of variables (H(t), p(t),p(t),G(t),A(t)). Many 
relations have been proposed in the literature from the purely phenomenological point of 



view [jH], 2C[. In the next section we discuss the possibility that the fifth missing equation 
is a relation of the form A = A(H). Most important, we will motivate this relation at 
a fundamental QFT level, specifically within the context of the renormalization group in 
curved space-time, and then discuss the possible form of the associated RG evolution law 
for the gravitational constant. 
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3. Renormalization group cosmologies 



Following the approach of J?], [2^, 23, 32] we describe the scaling evolution of A by in- 
troducing a renormalization group equation for the cosmological constant. The connec- 
tion between scaling evolution and time evolution is to be made in a second step. From 
dimensional analysis, and also from dynamical features discussed in the aforementioned 
references, the RGE may take in principle the generic form 

n 

Here \x is the energy scale associated to the RG running. Furthermore, in a RG cosmology 
with a time-evolving G we expect a RGE for G as well. This is best formulated in terms 
of G _1 as follows: 

^(I)^ B „„- (3.2, 

It will prove convenient for further reference to introduce the integrated form of these 
equations: 

A(/i) = ^ C n p 2n , = D n f n ■ (3-3) 

n " ' n 

Coefficients A n , B n , C n and D n in these formulas are obtained from a sum of the contribu- 



tions of fields of different masses Mj. The r.h.s. of (3.1) and (|3.2| ) define the (3- functions for 
A and G~ l respectively, which are functions of the masses and in general also of the ratios 
between the RG scale \i and the masses. For \i <C Mi, the series above are in principle 
assumed to be expansions in powers of the small quantities n/Mi |?], |2^]. However, the 
possibility to extend the analytical properties of some of these expansions will be discussed 
later. We assume that \x is of the order of some physical energy-momentum scale char- 
acteristic of the FLRW cosmological setting, and can in principle be specified in different 
ways. In the A-running framework of 0, it was assumed that p is given by the typical 
energy-momentum of the cosmological gravitons of the background FLRW metric. This 
suggests we can set // ~ H (Hubble parameter) , which is of the order of the square root of 
the curvature scalar in the FLRW metric, /j, ~ i? 1 / 2 . 

In the cosmological applications, e.g. when considering the evolution of the present 
universe, we are going to keep the Ansatz p ~ H. We will claim that it is a sound RG 
scale-setting also for a framework based on running A and G. The consistency of this 
choice will be discussed in Section || in the light of the scale-setting procedure devised in 
1 36]. For the moment we observe that, within the p ~ H Ansatz, the RGE ( |3.1|) above 
defines a relation of the type A = A(iT), and therefore provides the fifth equation suggested 
at the end of the previous section, needed to solve the cosmological model. At the same 



time Eq. (3.2) defines a relation of the sort G = G(H). Obviously these two functions of H 
cannot be independent, for the Bianchi identity imposes the differential constraint (|2.9| ) in 
which p = p(t) is a known function obtained after solving (|2.7l ) - once the equation of state 
P = p(p) f° r matter and radiation is specified. Put another way: consistency requires that 
the two expansions ([O]) and ( |3.2| ) ought to be interconnected via the differential constraint 
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2.9| ), so that if one of the expansions is given the other can be retrieved from the Bianchi 



identity. For example, in Ref. [23] we suggested the following specific form for the RGE 

'' /A ' a^ 2 M 2 + ..., (3.4) 



efln fj, (47r) 2 

with \i = H. Here we have defined the effective mass parameter 



M = 



1/2 



(3.5) 



and a = ±1 indicates the sign of the overall ^-function, depending on whether the 
fermions (a = —1) or bosons (a = +1) dominate at the highest energies. The form ( fO| ) 
for the RGE is based on the leading participation of the heaviest possible masses Mi in 
the expansion. These masses could arise e.g. from some Grand Unified Theory (GUT) 
whose energy scale Mx is just a few orders of magnitude below Mp. The form ([O]) results 
from truncating the series ( |3,1| ) at n = 1. The term n = 0, corresponding to contributions 
proportional to Mf, must be absent if we want to describe a successful phenomenology [23 



24]. Actually, from the RG point of view it is forbidden because in the physics of the modern 
Universe we probe a range of energies where \x is always below the superheavy masses M{. 
In addition, dimensional analysis and dynamical expectations tell us that the remaining 
terms are suppressed by higher powers of fi/Mi ^2], (23j. In the absence of experimental 
information, the numerical choice of a M 2 is model-dependent. For example, the fermion 



and boson contributions in (3.1) might cancel due to supersymmetry (SUSY) and the total 
/^-function becomes non-zero at lower energies due to SUSY breaking. Another option 
is to suppose some kind of string transition into QFT at the Planck scale. Then the 
heaviest particles would have masses comparable to the Planck mass Mp and represent the 
remnants, e.g., of the massive modes of a superstring. 

The obvious motivation for this RG scenario (2^, [}2| is that the r.h.s. of ( [3.4| ) is of 
the order of the present value of the CC. To assess this in more detail, let us integrate 
Eq. ([O]) straightforwardly. We find an evolution law for A which is quadratic in n and 
can be expressed, with the help of the notation fl3.3j), as follows: 



A(/x) = C + Ci fi 2 , (3.6) 

where 

C = A -^M 2 pf i 2 , C 1 = ^M 2 . (3.7) 

Within our Ansatz for the RG scale, fiQ stands for the present day value of the Hubble 
parameter, Hq, and the coefficient v is essentially the ratio between the squares of the 
effective GUT mass scale M and the Planck scale up to a = ±1, depending on whether 
bosonic or fermionic fields dominate in the mass spectrum below the Planck scale: 

a M 2 . . 

(3.8) 
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This parameter was first defined in Ref. ]^TJ, pj]. Irrespective of the precise meaning of v 
in the context of QFT in curved space-time, the latter reference suggested also how to use 
v to parametrize deviations of the evolution laws from the standard FLRW cosmology 5 . 
Equation ( |3.6[ ) is normalized such that A takes its present day value Ao for /i = fj,Q. The 
highest possible value of the heavy masses Mj is Mp, and therefore we typically expect 
that the effective mass ( |3.5[ ) is of the same order. However, due to the multiplicities of the 
heavy masses, we cannot exclude a priori that M > Mp. To be conservative we will take 

\u\ = i/ = — ~ 0.026 (3.9) 

as the typical value of v, corresponding to M = Mp. A very nice feature of this RG picture 
for A is that its evolution, being smooth and sufficiently slow in the present Universe, 
may still encapsulate non- negligible effects that could perhaps be measured 6 . In fact, the 
predicted rate of change of A (in other words, the value of its ^-function) at fj, = Ho is of 
order of A itself: 

0° A ~ C X H 2 \ t=t0 = |^ Ml H 2 < 10- 47 GeV* . (3.10) 

In contrast, at the earliest times when H ~ Mp, the same RGE predicts a CC of order 
A ~ Mp, which again is a natural expectation. If we define the energy scale E\(t) = A(t) ' 
associated to the vacuum energy density p\ = A at any given moment of the history of the 



Universe, the RGE above may explain [23, 32 1 why the millielectronvolt energy associated 



to the present value Ao happens to be the geometrical mean of the two most extreme energy 
scales in our Universe, Hq ~ 10 2 GeV (the value of \x at present) and Mp ~ 10 19 GeV: 

E A (t = t ) = Al /4 ~ VMpH = O (l(T 3 ) eV . (3.11) 

Amusingly this is also the order of magnitude of another, independent, cosmic energy 
scale: the one associated to matter density, EM(t) = /9^ 4 (t). Assuming k = and A = in 
Friedmann's equation (g~j), the energy scale associated to the matter density at present is 
seen to be Em (to) ~ ^JMp Hq. Remarkably the dark energy scale ( 3.11| ) associated to our 



RG cosmological framework coincides with the matter density scale defined by Friedmann's 
equation at t = to, and therefore it helps to explain the so-called "coincidence problem" [10], 
one of the two fundamental conundrums behind the cosmological constant problem. 

A value of v of order of 10~ 2 — 10~ 3 suggests that there might be some physics just 2 
orders of magnitude below the Planck scale, i.e. around 10 17 GeV. This would be indeed 
the mass scale of the heavy masses M, if we assume some multiplicity of order 100 in the 
number of the heaviest degrees of freedom in a typical GUT. To see this, rewrite (|3.8|) with 
v = 10~ 3 as 

d M i = 12 Ml v = 12 vr (1.2 x 10 19 GeV) 2 10" 3 ~ 5.4 x 10 36 GeV 2 . (3.12) 



5 See e.g. Eq. (4.11)-(4.14) of [g3| and discussions therein. Recently, this purely phenomenological point 
of view has been re-taken in ]42| where a similar parameter is called e (=3^). 

6 One could not achieve this extremely welcome pay-off if using RG scales other than fj, ~ H (e.g. 
/i ~ p 1,/4 ) and particle masses of order of those in the SM of particle physics. One unavoidably ends up 
with fine tuning of the parameters. See the detailed discussions in H 
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Then if we assume M{ = Mx, c% ~ 1 for all i, it follows that 

M x ~ 10 17 GeV for ^ c* ~ 100 . (3.13) 

i 

Therefore, there is the possibility that the origin of this RG cosmology could bare some 
relation to physics near the SUSY-GUT scale. 

From the previous considerations we realize that the models based on the RG evolu- 
tion of A could be serious competitors to e.g. quintessence models (and generalizations 
thereof) to account for some of the fundamental issues behind the multifarious cosmologi- 
cal constant problem ||, [TTJ] , and moreover they are also testable in the next generation of 
high-z SNe la experiments (such as SNAP |43]| ) - as has been demonstrated in great detail 
in The RG cosmological model analyzed in the last reference, although it constitutes 
a running A model, is a fixed-G model. In contrast, the RG model under consideration 
aims at preserving the virtues of the quadratic law ( |3.6|) for the A evolution (i.e. the rich 
phenomenology potential unraveled in ^4|) while allowing for a simultaneous evolution of 
the gravitational constant, G. This may add distinctive features and new phenomenologi- 
cal tests to be performed. Obviously it is important to analyze this possibility in order to 
distinguish between these two RG cosmological models and also to fully display the power 
of the RG approach versus the other dark energy models |H]] . 

4. On the RG scale-setting procedure in cosmology 

In this section we consider the consistency of the RG scale choice in certain problems related 
to the cosmological framework. Specifically, we discuss the consistency of our choice p, ~ H 
for the application to the large scale cosmological problems along the lines of the procedure 
recently proposed in [ 36 1 . In the next two sections we shall concentrate on the construction 



of a cosmological model based on the RG approach, and only in the last part of our work 
we will try to extend these methods to the local (galactic) domain. Let us commence with 
the scale setting at the cosmological level. Starting from the canonical constraint fl2.9| ) we 
trade the cosmic time for a generic scale parameter p (without specifying its relation to 
physics at this point), and rewrite it in the following way |36]| : 

Obviously, the matter-radiation energy density p is the only parameter which can not be 
defined as a function of p directly because it is unrelated to the RG in the first place. 
However, since G and A are primarily functions of p in the RG framework, the differential 
constraint (4.1) can be used to eventually define p as a function of p . If dp/dt ^ 0, it is 



£-0. (4,) 



easy to check that we can isolate p as follows: 

it \ 1 dA(p) ( d 1 V 1 fi . 

This scale-setting algorithm defines a relation p = f{p) between the cosmological data, 
represented by p, and the scaling functions of p. Under appropriate conditions, this relation 
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can be inverted to find /x as a function of the external data: /i = Since this function 

has been determined by the Bianchi identity, it should define a valid set of RG scales for 
cosmology. In other words, the specific scale /i to be used in cosmology should be chosen 
among the class of (potentially multivalued) functions obtained by this procedure under 
the condition that the scale so obtained is positive-definite and is a smooth function of the 
cosmological parameters. Further (physical) criteria may be needed to ultimately select 
one particular function within this class. 

Assuming that this is a valid method to define the RG scale, we are now in position 
to develop a scenario with simultaneous running of the parameters G and A in which the 
choice jj, ~ H is compatible with the aforementioned scale-setting procedure. To realize 
this scenario we stick to our RGE for A as given in ( |3.4| ) - hence to the quadratic running 
law ( |3.6|) - and we seek for a suitable RGE for G^ 1 . It cannot be an independent equation 
once the RGE for A has been established. Therefore, this restriction provides a method to 
hint at possible forms for the RGE ( |3.2|) compatible with the canonical constraint (|2.S|). 
For example we could start with an expansion akin to that in Eq. (3.1) of p3] for A: 

< 4 " 2 (g) = £ « M - + £ « M ° (m) 2 + £ < M - (£) 4 + - <«' 

Since the cosmological RG scale [i is assumed to be very small for the present universe, 
and the series is supposed to be convergent, we expect that the first term in this expansion 
is the crucial one and the remaining terms are inessential. Indeed, selecting the first term 
we are led to a RGE of the type 



w\a)-" M * + - (4 - 4) 

Consider now the generic expansions ([3.1|)-( |3~3] ) for A and G _1 . In particular notice the 
following relations among the coefficients for n = 0: 

C = A + A In ii - - A! $ - - A 2 $ + - 

D = + B In JL - \ £?x & ~ \ B 2 14 + ... , (4.5) 

where is the normalization point corresponding to the t = to moment of time. The 
coefficients (|4.5|) play a role in computing explicitly the r.h.s. of ([4.2|). In particular notice 



that Ao = 0, to be consistent with (3.7)(cf. Ref. p3[] ). Also important is that the scale 
derivative of Dq is non- vanishing, 

dfj, H 07T n 



Here we have used -Bo = 2z/Mj, = aM 2 /Qir as seen from (fO|) and ( |3.8[) . This result 
shows that G~ l is actually non-analytic at // = 0, a fact to which we shall return later. 
With these relations in mind, a straightforward calculation of the r.h.s. of Eq fl4.2|) , using 
the expansions ( |3.3| ), leads to 

2A) A 2 



p^-Co + dl-^-l) + - (4.7) 
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where we have neglected the higher order terms in /x. Next we substitute the coefficients 
( |3.7| ) in (|4.7[) and use the fact that Do — Mp. Irrespective of the sign a = ±1, the resulting 
equation determining /i is amazingly simple: 

1 - Cl (^-l) - 3 Mp ~ (48) 

where the last equality makes use of Friedmann's equation ( |2.4| ) in the flat case. Our 
starting Ansatz /i = H is, therefore, explicitly borne out within the criterion defined by 
Eq.fJ) (q.e.d.). 

The previous proof relies on the expansions (jO|)- (l3.3| ) advocated in the literature @, 



22| , |23| , pq ]. Notwithstanding, this does not exhaust the list of possibilities made available 
to us by the Bianchi identity in the canonical form fl2.9|) . It is still possible to extract 
the same phenomenology if we make allowance for an extended RGE for G beyond an 
expansion in (integer) powers of the small parameters /i/iWj. Then other (non-analytic) 
contributions may arise. We shall discuss this possibility along the second part of the next 
section. 

5. A cosmological model with logarithmic running of G and quadratic 
running of A 

In the previous section we have established a scenario where the choice fi ~ H is a consistent 
scale setting for cosmological considerations. Let us use this RG framework and apply 
it to a spatially flat Universe, i.e. k = in Eq. (|2,4|). This situation is not only the 
simplest one, but it corresponds (as already mentioned in the introduction) to the preferred 
scenario in the light of the present data from distant supernovae and the CMB. Moreover, 
the k = Universe already captures all the main traits of the new RG cosmology with 
variable G without introducing additional complications related to curvature. By collecting 
Friedmann's equation, the quadratic RG evolution law for A, and the (canonical) differential 
constraint imposed on the functions G(t) and A(t) by the Bianchi identity, we can set up 
the system: 

. 3H 2 
P + A-- 



lirG ' 

A = C + C 1 H 2 , (5.1) 
(p + A) dG + Gdk = 0. 

This system can be immediately solved for the function G = G(H; u) depending of the 
parameter u, with the result 

where G(Hq) = Gq = 1/Mp. In this way we have obtained explicitly the corresponding 
running law for G as a function of the RG scale [i = H. From this law we learn that for v > 
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the gravitational constant G decreases logarithmically with H, hence very slowly, which 
is an extremely welcome feature in view of the experimental restrictions on G variation 
(see later on). Such a decreasing behavior with the cosmological energy scale \x = H is 
equivalent to say that the law that we have obtained nicely corresponds to an asymptotically 
free regime for u > 0, much in the same way as the strong coupling constant in QCD. In 
addition, we discover that the parameter v defined in ( |3.8[ ) plays essentially the role of the 
(dimensionless) /3-function for G. Finally, we note that the function ( |5.2[ ) is non-analytic, 
which was expected from a remark made in Section ^. 

Consistency with the Bianchi identity requires that Eq. ( |5,2| ) must be the solution of 
the RGE (|3.2|) for G~ l . To check this, we note that in leading order the last equation boils 



down to Eq. (4.4). Integrating this equation under the boundary condition G(fi = /io) = Go 
we immediately recover Eq. ( |5.2|) , as desired. Had we ignored the first term in ( ^4.3[) and 
assumed that the first leading contribution comes from the second term (i.e. the one 
proportional to /i 2 = H 2 ) we would have obtained 

instead of ( |5.2| ), where £ is a dimensionless coefficient. The resulting scenario is not very 
appealing as it does not lead to any phenomenological consequence to speak of. In fact, 
if we substitute (^J^) in (|2.9| ) and integrate the resulting equation, the following evolution 
law for the CC ensues: 

A(H) =A + ^-(H*- F 4 ) . (5.4) 



Although this quartic law is - as the quadratic one (3.6)- smooth and mathematically pos- 
sible, it leads to absolutely no phenomenology due to the smallness of the Hubble parameter 
at the present time: Hq ~ 10~ 42 GeV. Similarly, the square dependence Go H 2 = H 2 /Mp 
in the running law for G given by (|5.3j ) is extremely weak. Therefore, it is clear that only 
if the first term of the r.h.s. of (|4.3| ) is non- vanishing it is possible to meet the logarithmic 
law ( |5.2[) for G and the quadratic law ( |3. 6| ) for A. This combination of RG laws is the only 
one compatible with the canonical constraint ( |2.9[) that may lead to a cosmological model 



with a rich phenomenology comparable to the model presented in |2^, 24]. 

A remark is now in order. As the reader could readily notice, by admitting that the 
first term of the r.h.s. of (|4.3|) is non- vanishing, we are departing from the decoupling 



theorem [45] and assume the non-standard form of decoupling for the parameter 1/G. In 



order to justify this departure from the canonical form of decoupling, let us remember 



that the attempt to derive the gravitational version of the decoupling theorem [37] was 
successful only in that sector of the gravitational action which is available in the linearized 
gravity framework. In the cases of A and 1 /G the explicit derivation of the decoupling law 
requires a new technique of calculation which should take the non-trivial and dynamical 
nature of the metric background into account. In the actual phenomenological approach 
it is primarily the Bianchi identity (rather than the decoupling theorem) that decides 
about the evolution of G once the evolution of A has been ordered. In this sense, the 
final structure for the RGE of G can be of the form (13.41) even if this does not fulfill the 
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standard expectations from the decoupling theorem. This uneven status between the two 
parameters can be understood on physical grounds, in the sense that the ultimate origin of 
G and A and their role in defining the background geometry may be essentially different. 
A can be viewed as a small parameter in QFT: this can be technically implemented by 
allowing for a fine-tuning of the renormalization condition for the vacuum counterpart of 
the cosmological term - see the Introduction. At the same time the situation with 1/G is 
quite different. For example, it could be that the Planck mass is not only some threshold 
scale for the quantum gravity effects but, after all, a real mass. One can imagine that the 
non-perturbative phenomena of quantum gravity are involved to generate this scale and 
then the standard decoupling law (which is essentially a perturbative result) would not 
apply to the parameter 1/G. 

We shall finish this section by presenting additional arguments concerning our choice 
for the identification of the scale parameter \x with the Hubble parameter H at low energies. 
These arguments do not look superfluous in view of the distinct choices for /i which one can 
meet in the existing literature plj , [?], 25, 36]. If one allows for some generalization 



of the expansions (|3.3|), it may have an impact on the scale setting. It is thus advisable 
not to assume a priori the result /j, = H derived in ( |4.8| ), although we want to show that 
under certain conditions we can also retrieve a similar scale setting. On physical grounds 
we cannot exclude a priori the presence of non-perturbative and/or non-local effects in the 



far IR [44]. Whereas at higher energies the behavior of G(^) -1 is influenced by the higher 
powers of \i in ( |3.3|) , the running in the present day Universe is controlled by the first 
coefficient Bq in (|3.2j), which in turn leads to a non-analytical behavior of Dq = -Do(^) 
as a function of \i in the series ( |3.3[ ). In other words, at the present IR energies the 
expansion G~ l in powers of \x may not be a pure Taylor expansion. Baring in mind these 
considerations, let us generalize the picture as follows. Assume that in the far IR the 
expansions ([T^) take contributions of the form 

A(/i) = a + ai (-^) , — !— = b + h ( ) , (5.5) 

where we assume arbitrary coefficients ao, b<j, ai,b\ and arbitrary exponents a, (3, which can 
be smaller than one. From these general expressions we may compute the corresponding 
RG scale following from (|4.2j). At the end of the calculation we encounter the result: 



ai 

P + a = j 



\HoJ h \fj, 



(5.6) 



This example already shows that even in relatively simple situations it may be difficult to 
find the explicit relation fi = f~ 1 (p) mentioned in Section |4|. It also shows that the scale 
setting fi ~ H will not always be possible. And, most important, it also helps us to find 
situations when it is actually possible. Indeed, let us consider a generic case where the 
RGE for G _1 takes terms of the form 
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For any 7 > (integer or not) this expression does respect the decoupling theorem, in 
contrast to ([4.3|). By integrating ( |5.7| ) we are led to an expression for G( / u) _1 as in (|5.5|), 
with 



Together with ( |5.7| ) let us assume that the A-scaling law is given by ( |5,5|) with a = 2. Then 
in the limit 7< 1 (hence (5 <C 1) we may solve approximately Eq. ( |5.6| ) for the RG scale, 
and we find 

2 2 /3bi p + a 8 ir rr2 . 

M - Mo ^— f—f - P + A = H 2 . 5.9 

2ai o +0i 3 Mp 

In deriving this formula we have used: i) the explicit form of the coefficients in Eq. (|5.S|); 
ii) ao = Co ~ Ao and a\ = ^qCi, which follow after comparing ([5l]) with ( |3.6| )-(|3?7D; iii) 
(3b\ ~ 2v Mp, which holds in the limit 7 <C 1. We have thus met again the scale /i ~ H. 
We see that it can be obtained within the more generic class of RGE's for G of the form 
( |5.7| ) in the limit of small 7 (corresponding to an effective logarithmic behavior), and when 
the cosmological term evolves quadratically with /j,, i.e. a = 2 in ( |5.5| ). This more general 
analysis may suggest the following interpretation. Recall that the cosmological domain 
has a natural finite size, which is defined by the largest cosmological distance where we 
can perform physical measurements and where all quantum fluctuations remain confined 
|38[| : to wit, the horizon dh ~ 1/H. The energy scale associated to the horizon is precisely 
our running scale \i = H. We do not consider graviton momenta below this scale. The 
underlying effective action may develop an IR tail of non-analytic contributions related 
to non-perturbative and/or non-local effects that cannot be described by a pure Taylor 
expansion. We have modeled this possibility with a RGE for G of the form ( |5 . 7j) with 
small 7. In general we expect non-analytic effects also for the local systems at, say, the 
galactic scale, which are always finite-sized by definition -see more details on this matter 
in Section |^. 

In summary, in this section we have shown that the quadratic law (|3.6| ) for A and the 
logarithmic law ( |5.2[ ) for G are a pair of RG scaling laws that can be consistently formulated 
in terms of the scale \i = H. This is tantamount to say that these running laws G = G(/x) 
and A = A(/i), when formulated in terms of = H, can be consistently transformed into 
time-evolving functions G = G{t) and A = A(t) that satisfy the original (time-dependent) 
canonical differential constraint ( |2.9| ) reflecting the Bianchi identity of the Einstein tensor. 
Furthermore we have seen that the RG scale for these laws can be derived from the recipe 
Q4.2j ) either by assuming the expansion ( |3.2| ) with Bo / as the leading coefficient, or from 
the exponential form ( |5.7| ) in the limit of small 7; in both cases it amounts to a logarithmic 
behavior of G = G(fx). As we have advanced in Section ||, a quadratic law for A is very 
important for potential phenomenological implications. In the following section we explore 
some of these implications for cosmology, together with the new feature introduced by the 
logarithmically varying gravitational constant. 
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6. Numerical analysis of the RG cosmological model 



The RG cosmological model under consideration is based on the system of equations ( |5.1| ) 
and the ordinary conservation law (|2.7| ) for matter and radiation. While the solution has 
already been exhibited in the form G = G(H;i>), A = A(H;i>) in Section || we want to 
present also the solution in terms of the redshift variable, z = (ao — a) /a, following a 
procedure similar to 



24 1 . This is specially useful for phenomenological applications, in 
particular to assess the possibility that this model can be tested in future cosmological ex- 
periments such as SNAP [43]. Therefore, we look for the one-parameter family of functions 
of the variable z at a given value of v; 



G 



G(z;u), A = A(z;v). (6.1) 

We are going to analyze first in detail the solution ( |6.lD for the present (matter-dominated) 
universe, where p = pm, PR — 0, and then make some comments on the implications for 
the nucleosynthesis epoch. 

In contradistinction to the cosmological model discussed in [^, the functions ( |6.1| ) 
cannot be determined explicitly in an exact form. However, an exact implicit function of 
the type F(G(z; v), pm(z); is) = can be derived after integrating the differential equations. 
For the matter-dominated epoch the result reads 



with 



1 + v In 



g{z\v) = 



1 - vg(z;v) 
,g{z;v){\ - v) 

Gn 



v In 



p M {z) + A - v p\ 



.01 



P°M + A 



Pm(z) = p° M (l + zy 



(6.2) 



(6.3) 



The value of the gravitational constant and the matter density at present {z = 0) are Go 
and p° M respectively, p° c being the critical density. We may check that ( |6.2| ) is satisfied by 
g(0; v) = 1 for any v, as it could be expected. The expression fl6.2j ) defines implicitly the 
function g = g(z; v), and once this is known the other cosmological function A = A(z; u) is 
obtained from 

A + vp M {z)g{z\v) 



k(z;v) 



vpl 



(6.4) 



l-vg(z;v) 

Obviously this satisfies A(0; v) = Ao for any u, as it could also be expected. 

Although the explicit analytic form cannot be obtained for arbitrary u, it is possible 
to derive an explicit result for small v -C 1 within perturbation theory. This situation 
actually adapts very well to the theoretical expectations on v based on effective field the- 
ory arguments 0, 23 1. However, there are also phenomenological reasons, mainly from 
nucleosynthesis [23, 24] and also from the CMB and LSS J|2|, leading to the same con- 
clusion 7 . Expanding the exact formulae above up to first order in v we may compute the 
relative correction to the gravitational constant: 



5g{z;v) ee 



G{z-v)-G 
Go 



1 ~ —v In 



p M {z) + A 
Pm + A o 



-v in (n M (z) + n° A ) 



(6.5) 

7 For instance, in [E3 it is shown that values of v (e/3 in their notation) of order of vq (Eq. (3.9)) or less 
are compatible with current data on CMB and LSS. 
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Figure 1: 



The dependence of the gravitational coupling constant G = G(z; v) on the cosmological 

We assume a 



redshift z for the six values (6.8) of the parameter v, where v§ is given by Eq. (3 
flat cosmology, k = 0, with = 0.3, fi^ = 0.7. For the parameter values of the order i/q, there is 
a significant variation of G at the redshifts amenable to high-z SNe la observations. 



where Qm{z) = pm(z)/ ' P° c - Recalling that the present day cosmological parameters, £l Q M = 
Pm I Pc an d = ^-o/Pc = Pa/Pc i do fulfill the sum rule Q° M + f2 A = 1 for the flat space 
geometry, we see that the previous expression satisfies Sg(0; v) = 0, as expected. The 
above formula can also be obtained from Eq. ( |5.2p at first order in v upon replacing H 2 
with the expression dictated by Friedmann's equation ( |2.4| ) with k = 0. For v > (resp. 
v < 0) the gravitational constant was smaller (resp. higher) in the past as compared to the 
present time. This was already patent from Eq. (|5.2[) when we remarked that, for v > 0, G 
becomes an asymptotically free coupling, which means that at higher energies (i.e. when 
we look back to early times) the coupling becomes smaller. Similarly, the explicit result 
for A(z; v) can be cast, to first order in v , as follows: 



5A(z; v) A(z; u) - A fi° 
For convenience we have expressed the result in terms of the relative deviation of A(z; v 



s ^=^jdr= I - ^ [(!+»)■ -i]. (6.6) 



with respect to the present value Ao, in order to better compare with Ref. 23, Amaz- 
ingly, to first order in v, the cubic law fl6.6|) is identical to the one obtained in the latter 
references where a RG cosmological model with variable A and fixed gravitational constant 
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Figure 2: The dependence of the cosmological term energy density A = p\(z; v) on the cosmo- 
logical redshift for values of the parameter v as in Fig. |]. For these parameter values there exists 
significant variation of A (expressed in units of its present value A = p A ) which may be detected 
by the high-z SNIa observations. 



was considered. This means that the promising phenomenological analysis made in [24] 
concerning the redshift dependence of the cosmological constant remains essentially intact. 
Let us, however, mention that the two models depart from one another to second order in 
v. For the model under consideration we find 



t|[(i + -) 3 



1 (u + i/ 



in [n° M (i + zf + ^ 



(6.7) 

To first order in v we recover equation ([O]), and moreover we see that <5a(0; v) = for any 
v, as it should be. 

The previous (approximate) results are confirmed and refined through a numerical 
analysis of the exact formulas ( |6.2| )-( p^4] ). We take as representative values of v some 
multiples of the typical value z^o, Eq. (|3.9|) . Specifically, we consider the six cases with 
different signs represented by the i^-data row: 



1 1 



(6.8) 



The results are presented in Fig. |l] and |2] for the standard values of the cosmological parame- 
ters £l° M = 0.3, 0^ = 0.7. Specifically, in Fig. || we plot the exact function G = G(z; v) 



for the six values of v listed above. Similarly, in Fig. |^ we plot the function A = A(z; v) for 
the same values of the v parameter. We have verified these results also by direct numerical 
solution of the differential equations using Mathematica [47]. By looking at Fig. || we may 
read off the relative correction 5q on G. We see, for example, that significant corrections 
5g — —5% are predicted at redshift z = 1.5, if v ~ vq. Furthermore, for the same values 
of v and z, we find <5a — 16% (cf. Fig. [2]). This correction to A should be perfectly mea- 
surable by SNAP [43]. Recall that this experiment aims at exploring cosmological depths 
up to z = 1.7, and that other experiments are planning to go even beyond (including the 
HST) @. 

We should point out that cosmological corrections to G of a few percent at high 
redshifts are not in contradiction to present bounds on the gravitational constant |49| , |50| . 
The reason is that the existing bounds are "local", i.e. at the scale of the Solar System, 
so that they do not provide any information on possible variations of G at cosmological 
distances. Actually it is difficult to measure the correction 5q at the cosmological level, 
at least in the present matter-dominated epoch. There is no direct way to measure G at 
distances of hundreds to thousands of Megaparsecs, unless we re-interpret the observed 
cluster dynamics in terms of G variation on a cosmological level, rather than in terms of 
CDM. We will further elaborate on this kind of idea in the next section, but restricting 
ourselves for the moment to the more familiar galactic domain. On the other hand, if we 
extrapolate back our running formula for G to the radiation epoch we might expect to find 
some additional information. Particularly relevant is the nucleosynthesis epoch (z ~ 10 9 ). 
Performing the corresponding change pu —> PR in Eq. (|6.3| ), with pr(z) = p° R (1 + z) 4 , 
we find that for v of order 10 -2 the relative variation 5g lies within the limits set by the 
latest analyses on G variation at nucleosynthesis, which give a margin of ~ 40% at the 95% 



C.L. [51] Thus the existing nucleosynthesis limits placed on v, based on considerations on 
A variation alone |23|, |1] , are essentially unchanged. 



7. A possible application to flat rotation curves of galaxies 

In this section we speculate on the possibility to apply similar RG methods to the more 
restricted astrophysical domains, such as the galactic scales. While we have seen that 
it is difficult to measure the correction 5q at the cosmological level, let us explore the 
possibility to measure local variations of the gravitational constant (induced by the RG) at 
the level of a typical galaxy, and extract some potential consequences concerning the origin 
of the flat rotation curves in its halo - usually associated to the presence of undetected 
CDM Q 8 . Here again the main question is to identify the scale of the RG. It can no longer 
be p o± H because measurements within our galaxy are not cosmological measurements, 
just astrophysical. Indeed, H does not change on the scale of a galaxy (a bound system!). 
However, if G varies according to some RGE, there must be some new order parameter for 

8 We point out that previous attempts exist in the literature to use the RG to explain the flat rotation 
curves of galaxies |32|. However, these old papers are based on the RGE of higher derivative quantum 
gravity. In this theory the RGE's for individual A and G _1 are gauge dependent and only the RGE for a 
dimensionless ratio is unambiguous - see e.g. 123] and references therein. 
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the local RG scale. For simplicity let us treat the galaxy as a central core surrounded by 
a spherically symmetric mass distribution. Then, from simple dimensional arguments we 
may assume that the appropriate RG scale for the galactic dynamics is given by 

/x = £, (7.1) 

where k is some (dimensionless) constant and r is the radial distance (in the "renormalized 
metric" , see below) measured from the center of the galaxy. If we assume that the renor- 
malization effects from the RG are small enough (after all they cannot change dramatically 
the Newtonian picture), we may take r as the standard radial coordinate in classical ce- 
lestial mechanics. Notice that p can never go smoothly to zero because r is bounded from 
above. This means that at the galactic level we do not expect that the corresponding 
expansion G(p)~ 1 in ( |3.3j ) is analytic near p = because the system is intrinsically finite- 
sized. However, by artificially extending the size of the local system we expect that the 
same IR behavior of the cosmological case mentioned in Section [| should be encountered. 
This provides a connection between the two regimes (local and cosmological) that can be 
used to establish the consistency of the RG application to the local case. The Ansatz 
( |7.1| ) obviously implies that spatial gradients of G (and may be also of A) are allowed. It 
should be clear that this is not in conflict with the Cosmological Principle because the 
latter applies to the Universe in the large, not to local systems such as galaxies or clusters 



of galaxies . Although the Ansatz (|7.l| ) is plausible, we may try to motivate it from the 



point of view of the scale setting procedure (4.2). Strictly speaking, the latter has been 



rigorously formulated only for the transition from the cosmic time evolution to scale evo- 
lution. However, given the fact that at very large scales r is of order of the horizon, 1/H, 
we see that with the hypothesis ( |7.1|) the original scale setting used in previous sections, 
p ~ H, is consistently met on cosmological scales. Consider the typical density profile for 
flat rotation curves: 

Here ro is a reference galactic length of the visible part of the galaxy; e.g. in the case of our 
galaxy it can be taken as the "local" position (at the galactic scale) of our Solar System, 
where the gravitational constant is supposed to be Go essentially everywhere, and po is the 
average density of the galaxy well before reaching the halo regime. Clearly p(r) — > po for 
r <C whilst p(r) falls off as 1/r 2 for r > ro (i.e. deep in the halo). This behavior is 
exactly what is needed to have a linear growth of the galactic mass deep in the halo region, 

M{r) = I d 3 r p(r) (x r , (7.3) 
and therefore a constant rotation velocity in it: 



GM(r) 

v = \ — -> const. (7.4) 



9 The Cosmological Principle indeed applies only to fundamental observers, and therefore to extremely 
large scale systems at the level of superclusters of galaxies pa] . 
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This is the classical "explanation" for the flat rotation curves, once the galactic density 
profile ( |7.2[ ) is accepted by fiat. One possibility (the most commonly accepted one) is that 
( |7.2| ) is associated to the presence of CDM distributed in the halo region r >ro. However, 
an alternative justification for this density function (and perhaps also the ultimate origin 
of the flat rotation curves) could actually be linked to the RG running of the gravitational 
constant as a function of the new scale parameter (7.1). A hint that this could be the 



case is obtained by substituting the density ( |7.2| ) into p in Eq. (^). It is easy to see that 
when the halo regime sets in, the resulting formula behaves precisely in the form ( |7. 1|) . We 
will see below that the density ( [7.2|) although it is not the real mass density of the galaxy 
in our framework (because we do not assume a priori the existence of CDM), it can be 
reinterpreted as the "renormalized" mass density, i.e. the "effective mass density" of the 
galaxy after considering the renormalization effects on Newton's gravitational constant. 
In this sense it is this expression that must be used to determine the RG scale. Let us 
therefore take these features as heuristic motivations to adopt [i from ( |7.1| ) as a valid RG 
scale for the galactic regime, and let us explore the ensuing consequences within our RG 
framework. Upon replacing /x = H with /i = n/r and Hq with k/tq in our running formula 
flgg ) we find 

1 + u In (?o/r z j 



We assume v > in (7.5) so that G is asymptotically free, and therefore it increases with 
r. This kind of growing behaviour of G with the distance is reminiscent of the above 
mentioned possibility that gravity can lead to non-perturbative effects in the far IR. The 
running of G(r) in the surroundings of the Solar System (i.e. for r ~ r$) is very small; in 
practice there is no running at all, for the entire Solar System is represented by the point 
r = ro without extension. This insures that we normalize our Gq above with the value 



obtained from highly accurate measurements of G performed in our neighborhood [|48|, 50]. 
Any departure from this value was never accessible to us because we never measured G 
(at least directly and accurately) outside the Solar System. Moreover, the formula above 
cannot be applied for r near r = 0. The reason is that the geometry is governed by the 
Schwarzschild metric (see below), and therefore there is an event horizon at some point 
r = 2GA4 <C ro, where M is the mass of the galaxy concentrated at its center. We 
can only assume valid Eq. ( |7.5[ ) well beyond this limit. Next we consider the modified 
Einstein-Hilbert action with variable G: 



Seh = ~ I d 4 xV=g. ia _^ ■ (7.6) 



B_ 

16?rG(r) 

If G were constant, we would expect that the metric that solves the galactic dynamics at 
large r is the spherically symmetric Schwarzschild metric. For variable G, however, the 
situation is more complicated. One could vary the action ( f7.6| ) and try to find out the field 
equations, but it is hard to solve them exactly. Whatever the solution is, at least we know 
that the spherical symmetry must still be there. Therefore, in the absence of an exact 
result we may use the fact that G(r) is given by ( [7.5[ ) and that v is a small parameter. 
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Following the conformal symmetry methods of [54, 13 1, we introduce the auxiliary action 

with respect to the auxiliary metric g^ u , and with constant G = Gq. This action transforms 
into the action 



o(3) 



±— [ d 4 x^j e 2s \R + 6{ds) 2 ] . (7.8) 
irG J 



' EH ~ 16vrG 
under the local conformal transformation 

9^ = e 2s g^v , (7.9) 
where s = s(r) is a function of r. With this trick, let us now choose the following conformal 

e 2s = = 1 + - In (r 2 /r 2 ) , (7.10) 

where G(r) is given by (|7.5| ). Then the action (|7.8|) becomes the r-dependent G action 
J7.6Q except for the extra derivative terms (ds) 2 : 



16vrG(r) 7 v 3 167rG(r) 

Notice that the square derivative terms are of higher order in v with respect to the La- 
grangian density yf^-g R/G(r), and can be neglected. Indeed, differentiating ( [7.1(30 on both 
sides gives e 2s ds/dr = —vjr and hence the square derivative that appears as the extra term 
in ( 7.11 ) becomes of order (ds) 2 ~ v 2 jr 2 . For a typical radius r comparable to the galac- 
tic size we can fully neglect this term as compared to R ~ 1/r 2 because contributions of 
order v 2 are very small. It follows that the solution of the field equations associated to the 
original action ( [7.6[ ) with variable G is, to first order in u, 

9^u = e~ 2s g^ u = (l + 2u In -0 g^ . (7.12) 

From the arguments above, g^y must be the solution of the ordinary Einstein equations 
with constant Go corresponding to the galaxy system treated as a point-like core of mass 
M. Therefore g^ on the r.h.s. of ( [7.12[) must necessarily be the standard Schwarzschild 
metric corresponding to constant Go: 

^-(l-^rf-j-^-^. (7,3) 

From ( |7.12|) we learn that the solution of the field equations for G = G(r) must be a 
modified (or "renormalized" ) Schwarzschild's metric. To first order in v the renormalization 
amounts to the new line element 

ds 2 = fl + 2Wn^ d4 chw . (7.14) 
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Let us now take the goo component of this modified metric and recall that the classic 
potential $ is identified from goo = 1 + 2$. To first order in v it reads 



M ( M\ r 

$( r ) = -G — + u 1-2G Q — In— . (7.15) 
r \ r ) ro 

The corresponding force per unit mass becomes 

d$ ~ x „ M fl „ M , r 2 



F(r) = -— = -(1 - 2 „) Go w - v I - + G ^ In ^ ) . (7.16) 



This intriguing result shows that for sufficiently large r (namely, in the halo of the galaxy 
where r >> ro) the force that we have found is not purely Newtonian as it does no 
longer decay as 1/r 2 at high distances, but as 1/r. The coefficient of the 1/r term is 
precisely the parameter v of our renormalization group model. Equating this force to the 
centripetal force —v 2 /r and neglecting for the moment the 1/r 2 (Newton's) term and the 
(In r)/r 2 term, one finds that for large r the rotation velocity is virtually constant, i.e. it 
does not depend on r (the rotation curve becomes flat). Such asymptotic velocity is thus 
determined in good approximation by the square root of our original parameter v = \fu. 
In conventional units the velocity reads v/c= \fu . Of course we expect v/c « 1 for any 
rotating object in the halo. But this can be well accomplished through the coefficient v, 
which was always expected small in our effective field theory approach to the /^-function 
of A down the Planck scale, see (|3.8|) . We argued that v could typically be of order 10~ 2 
for the cosmological considerations. However, we may ask now: how much small must it 
be to describe the physics of the flat rotation curves of galaxies? The rotation velocities in 
most galaxies are of a few hundred Km/s, say 200 — 300 Km/s Q. It follows that if this 
model is to describe them, v should be of order 

v ~ 10~ 6 . (7.17) 

Obviously a very small value! It implies that all of the approximations that we have made 
for v <C 1 should be perfectly acceptable. And of course, from the cosmological point of 
view, it is neatly compatible with the CMB and LSS measurements [^|. With this value of 



v we may explicitly check that the dominance of the leading extra term v jr in ( 7.16j ) over 
the Newtonian one, Go M/r 2 , really takes place in the halo of the galaxy, where r >> ro- 
Indeed, for a typical galaxy with a mass 10 10 " 11 times the mass M & of our Sun, and with 



the v value (7.17) obtained in our RG framework, it is easy to check that the v jr term is 
comparable to the Newtonian force just around a scale of several to ten Kpc, namely near 
the onset of the halo region. Beyond 10 Kpc, and of course deep in the halo region, the 
V- force" is dominant. 

If we keep the other terms in ( |7.16| ), except for the very small correction to the New- 
tonian part of the force, we find the following formula for the flat rotation velocity: 

v/c=a u + Go— (i + v ]n^) ~Jv+G — , (7.18) 
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where in view of the small value ( |7.17| ) we have neglected v ln(r 2 /rg) <C 1. The previous 
result reads essentially v/c = \fv up to a small correction which depends on the total 
mass of the galaxy and the radial distance to the rotating object in the halo. This is 
again an interesting feature, because the velocity of the flat rotation curves is well-known 
to show some departure from universality jl]] . In fact, we can further shape this formula 
such that it becomes closer to the observed behavior. Prom the modified Newton's law we 
may estimate the effective mass (or "renormalized mass «M r ") of the galaxy associated to 
the renormalization of G. By equating the purely Newtonian term in ( 7.16| ) to the leading 
correction proportional to v we immediately conclude that the "effective additional mass" 
of the galaxy is proportional to r in the halo region, 

5M(r;v) ~ (j^J r (r»r ). (7.19) 

The total renormalized mass of the galaxy (the one which is actually observed, if one 
adopts the purely Newtonian approach) is thus given by the "bare mass" (i.e. the "real" or 
ordinary mass A4 made out of baryons) plus the renormalization effect computed above: 
A4 r = M + 5M. It should, however, be clear that this is a huge renormalization effect, 
because 5M is comparable to (actually larger than) M as it could be expected from the 
fact that, in the language of the CDM, the extra amount of matter in a typical galaxy 
is known to be dominant over the luminous one: roughly ten times the total baryonic 
matter (whether luminous or not). It means that we should have 5M ^> A4, equivalently 
M r ~ 5A4. This feature is nicely realized in our RG framework and again emphasizes the 
potential non-perturbative character of the gravitational effects in the IR regime, similar 
to the strong low-energy dynamics of QCD. To see this quantitatively, let us take r deep in 
the halo region, say in the range 10 — 50 Kpc, and let us assume a galaxy with a physical 
mass A4 = 10 10 M & (typically the case of the total baryonic mass of our galaxy). It is easy 
to estimate from ( |7.19| ), using the value of v given in ( 7.17| ), that 8Ai ~ 10 Ai, which is 
in the right order of magnitude to describe a total effective mass around Ai r ~ 10 11 M & 
without invoking the existence of CDM. Strictly speaking, there is no physical halo in the 
RG approach, but at some point the local RG description must break down. This should 
define an "effective RG halo" of the galaxy, beyond which a larger scale RG picture should 
take over. 

The linear growth of the mass given in Eq. ( |T. 19|) implies that the "effective density" 
of the galaxy in the halo region behaves as 

P( r ) ~ (T=7r) (zs) ( r » r o) • ( 7 - 2 °) 



This effective density replaces the material effects associated to the CDM matter density 
Q7.2D . Comparing both densities in the region r > ^ we obtain Tq po = v/AttGq. Setting 
Po ~ M/tq we have M/tq ~ v/Gq, which is indeed satisfied within order of magnitude 
for a similar set of inputs as before. Finally, since the renormalized mass of the galaxy 
is proportional to the galactic radius in the halo region, we conclude that the rotation 
velocity ( 7.18| ) can be cast as 

v(r;M)/c~ ^ [1 + e{r; M)} 1/2 (r > r ) , (7.21) 
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where e(r; A4) is a small correction that exhibits the departure of the flat rotation curves 
from universality. It reads 

e(nM) = 1 ^- ) , (7.22) 



with Ai v (r) = 6Ai(r; v) given by Eq. ( [7.19 ). Notice that A4 u (r) does not depend on the 



mass of the galaxy, but only on v and on the radial position of the rotating object. We have 
seen above that typically Ai u (r) ~ 10 M. for r ro. Therefore, we conclude that the flat 
rotation velocity becomes a function v = v(r;A4). In other words, the small corrections 
to universality depend on the mass of the galaxy and of the radial location of the object, 
another observed feature. Only further refinements of the RG approach and the galaxy 
model may perhaps shed some light on the detailed numerical predictions concerning the 
value of the rotation velocity and its dispersion. Our formula above can easily accommodate 
velocity dispersions of order 10%. However, we cannot expect perfection. Recall that we 
modeled the galaxy as a spherical distribution of matter while it is roughly disk-shaped, 
thus other geometric and dimensionful variables might enter (including the mass M of 
the galaxy) to refine our scale Ansatz ( [7. 1| ) . This could have some further impact in the 
departure of the flat rotation velocity from universality. 

We should recall that the flat rotation curves of the galaxies can be fitted within 
the Modified Newtonian Dynamics (MOND) pB| . This phenomenological model requires 
gravitation to depart from the standard Newtonian theory in the extra-galactic regime 
(r S> ro). The MOND potential takes the form 

M 



$( r ) = -G — + y/a GM In — . (7.23) 
r ro 

The structure of the extra-galactic modification (the second term) is such that one can 
fit both the flat rotation curves and the departure of these curves from strict universality. 
The first feature is achieved from the logarithmic form of the correction, and the second 
feature is implemented by allowing for the coefficient of the logarithm to depend on the 
product GA4, where M is the mass of the galaxy. Obviously the overall coefficient must 
be dimensionless, and therefore one is forced to introduce a new "fundamental accelera- 
tion" parameter ao in Eq. ( [7.23 ), which can be fitted observationally to a value of order 



10~ 8 cm sec~ 2 . The combination of these two features is all what is needed to simultane- 
ously account for the leading and subleading observational effects of the galaxy rotation 
curves without invoking the existence of CDM haloes: to wit, the flatness of these rotation 
curves and the so-called Tuny-Fisher's law, v 4 oc C, where v is the asymptotic value of the 



rotation velocity and C is the absolute luminosity of the galaxy [56]. This law expresses 
some departure of the flat rotation curves from universality and can be rephrased in terms 
of the mass of the galaxy. If one makes the standard assumption that the mass-to-light 
ratio, M/C, is essentially constant for spiral galaxies, then v 4 oc A4. MOND indeed imple- 
ments this relation in the form u 4 = ao G M. - as it trivially follows from ( |7.23j ). However, 
it must be clearly stated that in spite of the great phenomenological success of MOND, so 
far none of the relativistic gravitational theories proposed up to date (aiming to underpin 
that phenomenological model from first principles) can be considered as free of theoretical 
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and experimental inconsistencies, including the recent attempt in |]57| . In this sense, it is 
a welcome feature that we can derive the leading effect (the existence of the flat rotation 
curves themselves) from our RG framework. This is obvious when we compare our modified 
Newtonian potential (7.15) and the MOND potential ( |7.23| ), The two potentials have in 
common the necessary logarithmic term to account for the leading correction to Kepler's 
law entailing the flat rotation feature. However, the RG-inspired potential ( [7.15 ) does not 
lead to the (subleading) Tully-Fisher's effect, while the MOND potential contains the piece 
V^o G M. just ordered ad hoc to reproduce this additional feature. It could be that the 
RG-potential misses this piece due to the simplifications we have made to reach the leading 
result. We hope that further investigations will help to clarify this issue. In the meanwhile 
Eq. ( 7.22 ) above already suggests that some mass-dependent departure from universality 
is predicted in the simplest RG approach. 



With v as small as ( 7.17 ) the increase of G(r) from (|7.5| ) remains negligible, even 
at cosmological distances. Therefore, it is conceivable that G(r) essentially behaves as a 
constant, with no further significant variation up to the horizon. At this point we recall 
that the motivation for ( |7.5|) was, in contrast to the more formal one of ( |5.2|) , essentially 
heuristic. However, since at very large scales r is of order of the horizon, 1/H, the function 
G = G(r) merges asymptotically with G = G(H). In this way the original scale setting 
[i = H is retrieved on cosmological scales and there is a matching of the two scale settings 
in the asymptotic regime. Notwithstanding, we should clarify that we do not expect a 
sharp transition from one RG regime to the other; we rather expect the existence of some 
interpolating RG description that connects scales of order 100 Kpc (which characterizes 
the maximum expected size of the galactic haloes) with scales of order 100 Mpc (above 
which the large scale cosmological picture takes over). In the crossover region the effective 
density ( 7.20| ) should fall off faster than 1/r 2 . Unfortunately, it is not possible to be more 
quantitative here on the specific form of this transition, but we do expect it should be 
there. After going through this transition the running formula ( |7.5|) becomes ( |5.2| ) and 
the common, and essentially constant, value of G in the two scaling laws represents the 
cosmological value of the gravitational constant. The only measurable effect left in this 
scenario would be the aforementioned correction to Newtonian celestial mechanics, which 
could have an impact on dark matter issues at the galactic and even at larger (cluster) 
scales. But we are not going to push further this possibility here nor to contemplate the 
possibility of other refinements. In particular, in the above considerations we have ignored 
the potential running effects from the cosmological constant at the galactic scale, which 
were assumed to be small. We close this section by noticing that with v of the order of 
( 7.17 ), the GUT scale potentially associated to these quantum effects would be Mx ~ 
10 16 GeV rather than 10 17 GeV (see Eq. ( |3.12| )), therefore closer to the lowest expected 
SUSY-GUT scale. 



8. Discussion and conclusions 



In this paper we have further elaborated on the formulation of the renormalization group 
cosmologies along the lines of 0, 21, 22, ^4|, 32, 33, 36]. The main aim of these RG 
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cosmologies is to provide a fundamental answer, within QFT, to the meaning of the most 
basic cosmological parameters such as Hq, G, p° M and Ao = p\, and at the same time to 
try to understand their presently measured values and possible interrelations among them. 
In doing so we have found that the classical FLRW cosmologies become "renormalization 
group improved", and therefore they become slightly modified by quantum corrections 
associated to renormalization group equations that govern the evolution of the parameters 
A and G. This evolution behaviour, primarily a scaling one within the framework of the 
RG, can be inter-converted into a temporal evolution (i.e. in terms of the cosmic time or, if 
desired, in terms of the redshift), after an appropriate identification of the scale parameter 
[i of the RG. We have found that [/, ~ H(t) is a consistent identification for the class of 
RG cosmologies in which A evolves quadratically, and G evolves logarithmically, with (i. 
And we have shown that this combination of running laws is a most promising one from 
the point of view of phenomenology. The scaling (and redshift) variation of these two laws 
is very simple and it is controlled by a single parameter, v - cf. Eq. ( |3.8[ ). For v > 
the running law for the gravitational constant exhibits (logarithmic) asymptotic freedom, 
much as the strong coupling constant in QCD. For v <C 1 the quadratic evolution law 
for A turns out to be very similar to the one previously found in |^, |24| and therefore 
it can be tested in the future experiments by SNAP ^jj, provided v is not much smaller 
than vq ~ 10~ 2 . This value is presently consistent with data from high-z SNe la, CMB 
and LSS. Furthermore, from the point of view of the cosmological constant problem, we 
have seen that this RG framework can shed some light to the coincidence cosmological 
problem (p9 ~ p a M ?) because it provides a natural explanation of the puzzle of why the 
energy scale associated to the cosmological constant is at present the geometric mean of 
the Hubble parameter and the Planck mass: E\(to) ~ \JMp Hq. 

It is important to remark that the quadratic law ( |3.6| ) associated to our RG framework 
is not just a proportionality law A ~ H 2 , but a variation law: 5 A ~ H 2 . In contradistinction 
to the former, the latter actually introduces an additive term. This is because it originally 
comes from a RGE of the form dA/dlnH oc H 2 , which has to be integrated, see Eq. (|3.4|). 
Therefore the ensuing cosmologies are very different from the A ~ H 2 ones. Only the latter 
have been amply discussed in the past in many places of the literature |y| [2(J, including 



some recently raised criticisms 58]. In contrast the 5 A ~ H 2 law is, to the best of our 



knowledge, relatively new in the literature @, and its motivation within QFT gives it a 
credit at a much more fundamental level. Let us mention that some variation of this law 
has been proposed very recently (on pure phenomenological grounds) in [jffil , in the form 
dA/da oc dR/da, where R is the curvature scalar. However, since p = R 1 ! 2 is, as mentioned 
above (and in 0), a covariant generalization of p = H, it is essentially the same law and so 
the conclusions are similar. Therefore the RG cosmologies associated to 5 A ~ H 2 can also 
be formulated covariantly, without committing to the particular FLRW metric. But from 
the practical point of view (and also to better assess its physical meaning and potential 
applications [24, |33|]) it is better to formulate this law making use of the Hubble parameter, 
instead of curvature. 

The last part of our paper is more speculative. We have tantalized the possibility of 
extending the RG formalism to local cosmological domains such as the galactic level. After 
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arguing that the local RG scale could naturally be /i ~ 1/r, we have solved the modified 
Schwarzschild metric in the limit v <C 1 and derived a RG correction to Kepler's third law 
of celestial mechanics. We have hypothesized that this quantum correction could furnish 
an explanation for the flat rotation curves of the galaxies, provided v ~ 1CT 6 . Obviously, 
this possibility would exclude measurable effects on the cosmology side, where a larger v 
is needed, but it might hit square on one of the biggest astrophysical conundrums of the 
last thirty years. This result is curious enough and rather intriguing, and we have found 
interesting to include it as an additional potential application of our RG framework. All 
in all, we deem this last part of our work more speculative than the rest, because we have 
had to descend from the cosmological level to the (more complicated) local domain where 
cosmology is entangled with astrophysics, and other factors may come into play. Thus, 
more detailed theoretical considerations are needed to formally establish the RG effects at 
the galactic level, which could be modified by other factors such as the modeling of the 
galaxy and the details of the RG framework. For this reason, even if v is of order of vq, 



and therefore much larger than (7.17), the corresponding RG predictions at the galactic 
level are not necessary in contradiction with the experimental situation. The actual effects 
could be much smaller due to, say, the need to introduce some significant modification in 
the Ansatz ( |7.l|) , as noted in the last section. Finally, even if accepting that v is small 
enough for the RG to be able to provide some explanation of the flat rotation curves, it 
does not preclude the possibility that galactic dark matter is present in the haloes of the 
galaxies. It only shows that other, even more subtle, effects could creep in. Whether or not 
these effects are directly related to the RG must be further investigated, but already at the 
pure phenomenological level the possibility that a logarithmic law of the type ( |7.5| ) could 
be behind the flat rotation curves looks rather enticing. Obviously, even if it is not the 
sole explanation, it may help to reduce the total amount of dark matter needed to account 
for the flat rotation curves and also to explain some features related to the deviation from 
universality of the flat velocities. In any case it is clear that further experimental work will 
also be necessary to decide on the ultimate origin of the rotation curves of the galaxies p9|| . 
In general, even if letting aside the possibility of these local effects, our cosmological RG 
framework provides a sound alternative to explain the nature of the cosmological parame- 
ters and their potential connection to quantum field theory. It also offers a link between 
the low-energy physics characteristic of large scale Cosmology with the high-energy inter- 
actions of Particle Physics. Last, but not least, it constitutes a ready testable framework 
for the next generation of cosmological experiments. 

Note Added: Shortly after the e-print Archive version of this paper was submitted 
as hep-ph/0410095, some other works were submitted to the Archive - see M, M- which 



overlap significantly with the content of Section 7. These works also use the RG to study 
the flat rotation curves of the galaxies and arrive to conclusions very similar to the ones 
presented here. Since, however, their approach is different from ours we think it can be 
useful for the reader if we comment on the similarities and also the main differences between 
the two methods. 

First of all we recall that two Renormalization Group formulations for gravity are 
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present in the literature, namely the perturbative RG for curved space-time and the Wilso 
nian non-perturbative RG for quantum gravity. The present work is based on the first of 
these RG approaches (see ]^7|, for an introduction and references). The object of our 
interest is the quantum theory of matter fields on classical curved background. It is very 
important that this theory is renormalizable and therefore application of the perturbative 
RG to this case is a consistent procedure. Unfortunately, the consistency does not mean 
that this approach does not meet real difficulties. The existing well-established formalism 
is essentially restricted by the MS-scheme of renormalization. Within this approach one 
can learn how the observables dependent on the dimensional parameter fx, but the physical 
interpretation of this parameter always requires a special additional effort. Furthermore, 
the MS-scheme is reliable only in the UV limit and is not appropriate for the investigation 
of the relatively weak IR scaling, that is for the study of decoupling. An alternative mass- 
dependent scheme has been recently applied in curved space-time, but only in the particular 
case of linearized gravity |T?j] . In this case we can observe the Appelquist&Carazzone de- 
coupling theorem [45] to hold in the higher derivative sectors of the gravitational action. 
The present and previous |7], 22, 23, 24] works are based on the phenomenological assump- 
tion that a similar kind of decoupling behaviour takes place for the cosmological constant. 
One can find an extended discussion of this issue in the Appendix A of |24|] . 



On the other hand, the alternative works f6C| (exactly as the previous ones j25|, g6] , 
see further references therein) deal with the IR effects of quantum gravity. It is well 
known that quantum gravity is inconsistent within the perturbative approach and hence 
the appropriate version of the RG is necessarily non-perturbative. Ref. |6(| is based on 
a simplified version of the Wilsonian approach to this non-perturbative RG framework. 
The objective is the derivation of the functional of the flow effective action, depending 
on the cut-off parameter k. This parameter can be viewed as a direct analog of the 
parameter [i for the perturbative RG in the MS-scheme. The main similarity between 
the two parameters is that one also needs to properly identify k with certain physical 
quantity (e.g. the energy of the gravitational quanta or the inverse of the cosmic time 
t). Without such identification the theory has no physical output. The choice of the 
scaling parameter in both RG approaches is not constrained by the RG formalism itself, 
and in certain (sufficiently simple) situations it can be guessed on the basis of physical (e.g. 
dimensional) arguments. This is the case when the system has just one relevant physical 
scale. For instance, in the astrophysical case under consideration in Section [?] the galaxy 
is assumed to have spherical symmetry and then the single relevant physical scale is the 
radial distance r from its center. This leads to the successful scale choice (7.1). On the 
other hand, in the cosmological context, a possible way of constraining the scale setting 
is the procedure applied in Section |2[ This leads immediately to our choice fj, = H - cf. 
Eq. ([D|) - roughly consistent with the scale setting k ~ 1/t advocated in the alternative 
RG approach of | |60| |. However, even this regular procedure does not obviously provide 
universal solution for all physical situations (for example, in multiple scale systems or 
in the transition from the cosmological to the astrophysical domain). Therefore, the link 
between the RG running and the r-dependence in the astrophysical setting is an essentially 
phenomenological assumption for both RG formulations, and needs further investigation. 
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As we already mentioned above, in our RG approach we assume the standard form 
of decoupling for the cosmological constant. Similarly, the alternative RG approach - the 
Wilsonian non-perturbative formulation of the RG for quantum gravity - requires also 
of strong phenomenological input. In the Wilsonian version of the RG used in |K| the 
effective action depends on infinitely many couplings associated to powers of curvature and 
non-local invariants. Solving the exact functional RG flow equation is virtually impossible 
at present and in practice one truncates this effective action into a linear combination 
of powers of curvature, the simplest possibility being the powers one and zero, i.e. the 
Einstein-Hilbert term and the cosmological term. The coefficients of these terms define the 
flow trajectory (G(k),A(k)) projected onto this simplified finite-dimensional space. 

However, by virtue of this simplification the flow trajectory G = G{k) , A = A(k) 
cannot be derived from the RG flow equation, and therefore an Ansatz on the trajectory 
is to be made in general. For the galactic problem at hand, one typically dismisses the 
role played by the cosmological term and is left with proposing an (ad hoc) modeling for 
the G flow. The model trajectory explored in [ SO | reads G(k) = l/k q (with q an arbitrary 
parameter). In order to describe the flat rotation curves the authors of |5(]] have to take for 
q a very small value of order 10 -6 . In contrast, in a previous paper |26| it was assumed the 
particular value q = 2 in order to describe a RG cosmology with a (hypothetical) IR fixed 
point. Obviously this trial and error strategy shows that the state of the art in the RG 
flow effective action approach is not good enough at present as to allow for determining 
the flow trajectory G = G(k) in a fully dynamical way. It is nevertheless remarkable that 
after accepting the value 10 -6 the previous power law trajectory becomes effectively a 
logarithmic running law, which is entirely similar to the one we have obtained in Eq. ( |7.5| ) 
on a very different basis. 

Comparing the merits of the two RG approaches we can indicate two advantages (from 
our point of view) of the one considered in the present paper. First, we are dealing with the 
quantum effects of matter fields, which definitely take place. The necessary phenomeno- 
logical implication concerns only the technical (albeit serious) difficulty of evaluating the 
vacuum effective action of massive fields on a curved background. In contrast, the rele- 
vance of the quantum gravity effects and the possibility of their description without solving 
the fundamental problem of quantum gravity looks rather uncertain. This concerns both 
the method of Ref. |k| and the one of the old works (5^] where the perturbative higher 
derivative quantum gravity effects were considered (and similar results were obtained). 
Second, the present approach has at least the virtue of introducing a parameter, u, whose 
definition (namely, v is the ratio of the masses squared of the heavy matter fields versus 
the Planck mass squared, see Eq. (|3.§|)), naturally fits in the right order of magnitude to 
explain the effect. At then end of the day, since both approaches are essentially based on 
phenomenological input, the preference of one of them remains a matter of taste. However, 
the net outcome is that the running law ( |7.5| ) and the one found in are identical if one 
identifies q = 2 v in the limit of very small value of these parameters (the only relevant 
limit for the galactic dynamics). 

Let us also mention a recent modification of the RG approach, the so-called metric- 
skew-tensor-gravity theory proposed in |61]. Here a skew-symmetric field tensor F^xa is 
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introduced and coupled to matter. This model introduces a new particle, but the main 
contribution to the mass of the galaxy is not from the skew field energy density, but from 
the RG flow running of G (induced by the coupling of this field to ordinary matter), which 
leads to a modified Newtonian acceleration law that can be successfully fitted to the galaxy 
rotation curves. Again this is achieved for special RG trajectories for which gravity becomes 
a (slowly varying) "confining force". In this case the trajectory is again the (/-dependent 
one mentioned above, but for the original value q = 2. The price of returning to this value, 
which is associated to a (hypothetical) IR fixed point, is to postulate the existence of a 
new particle in the galaxy halo. The physical consequence is again similar to our running 
law for G, Eq. (f7.5|), which entails anti-screening for v > 0, much in the same way as the 
QCD force at small distances. 

We find suggestive that all these RG approaches, even though they differ in the par- 
ticular details, seem to converge to the same basic idea that the flat rotation curves of the 
galaxies could perhaps be understood without using the hypothesis of CDM. At the mo- 
ment these RG frameworks cannot account for subleading effects, such as the Tully-Fisher 
law. However, they may constitute a promising framework to eventually construct the 
theory of the flat rotation curves of the galaxies from the first principles of Quantum Field 
Theory. The necessary modifications to account for subleading effects can be a matter 
of improving the technicalities of the RG framework and of the scale-setting procedure, 
including a more accurate modeling of the galaxy. 
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